CERN-TH/98-85 
FERMILAB-PUB-98 /093-T 
IC/98/25 
Eep-ph/98U4253| 



Extracting CKM Phases and B s —B s Mixing Parameters 
from Angular Distributions of Non-Leptonic B Decays 



Amol S. Dighe 

The Abdus Salam International Centre for Theoretical Physics 

34100 Trieste, Italy 

Isard Dunietz 

Theoretical Physics Division, Fermi National Accelerator Laboratory 

Batavia, IL 60510, USA 

Robert Fleischer 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



Abstract 

Suggestions for efficiently determining the lifetimes and mass difference of the 
light and heavy B s mesons (Bg,B^) from B s — > J/tpcj), D* + D*~ decays are given. 
Using appropriate weighting functions for the angular distributions of the decay 
products (moment analysis), one can extract (Th,^l, Am)B s . Such a moment 
analysis allows the determination of the relative magnitudes and phases of the 
CP-odd and CP-even decay amplitudes. Efficient determinations of CP-violating 
effects occuring in B s — ► J /il)(j),D* s + D*~ are discussed in the light of a possible 
width difference (Ar)s a , and the utility of this method for B — ► J/ipK* ,D* + D 
decays is noted. Since our approach is very general, it can in principle be applied 
to all kinds of angular distributions and allows the determination of all relevant 
observables, including fundamental CKM (Cabibbo-Kobayashi-Maskawa) param- 
eters, as well as tests of various aspects of the factorization hypothesis. Explicit 
angular distributions and weighting functions are given, and the general method 
that can be used for any angular distribution is indicated. 
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I Introduction 



Strategies for obtaining experimental insights into CP violation and non-factorizable con- 
tributions to weak decays are of particular interest for present particle physics. The ob- 
servables of angular distributions can be obtained in an efficient way by using an angular 
moment analysis ||J]— 0. In this approach, the observed experimental data are weighted 
by judiciously chosen functions, which project out any desired observable. This strategy 
is an alternative to the usual likelihood fit method M. It is demonstrated that the mo- 
ment analysis extracts all observables of measured angular distributions, such as the ones 
occurring in weak decays of pseudoscalars [P — > V£v, Xjiu, VV, etc.]. This method is 
of general validity. In our present paper, we apply the formalism to angular distributions 
U of B s and B meson decays into two vector-meson final states that are caused by 
b — > sec quark-level transitions. By making use of the general formalism outlined in this 
paper, it is straightforward to derive weighting functions for other exclusive mesonic or 
baryonic transitions, governed for instance by b — > cud, civ, u£v, c — ► sdu, s£ + u, d£ + v. 

The mixing between neutral B s mesons is expected to give rise to CP-even (5 S L ) and 
CP-odd (Bf) mass eigenstates, which may have a perceptible width difference Ar = 
T H — T L 0. Using appropriate weighting functions for the angular distributions of the 
decay products in the transitions B s — > J/ipcf) and/or B s — > D* + D*~, one can extract 
(T H ,T L , Am) Bs . 

A characteristic feature of the angular distributions considered in this paper is the 
fact that they contain terms describing interference effects between CP-even and CP-odd 
final-state configurations. Because of the lifetime difference, these contributions give rise 
to a term in the time evolution of the untagged rate, which is proportional to || : 

^Ht _ D -T L t 



(e-^-e-^jsin^cKM, (1) 

where 0ckm is a weak phase that is introduced through the CKM matrix ||. In the B s 
decays considered in the present paper, 0ckm is related to the Wolfenstein parameter 
f] |10|| . It is a remarkable feature that time-evolved untagged data samples of angular 



distributions of B s decays may exhibit CP- violating effects, if Ar is sizeable || |TTJ . This 



feature may be important, because it provides an alternative to previous investigations, 
which have shown how to extract sin </>ckm from tagged, time-dependent analyses |L2], p~3|j . 
This extraction, however, may not be feasible in the near future because it requires 
tagging and superb vertex detectors, which must resolve the rapid Ami oscillations. In 
contrast, any dependence on Amt cancels in untagged data samples, which therefore 
allow feasibility studies with current vertex technology \Ti\. 



Concerning tests of the factorization hypothesis |15[— pO], we divide the b — > sec 
modes into the following two categories: 

• colour-suppressed decays: B q — ► J/tpV with (q,V) G {(s, 4>)\ (d, K*°); (u, K* + )} 

mm- 



colour-allowed decays: B q — > D* + D* with q e {s, d, u] [|3|, |24 . 
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Whereas the validity of the factorization assumption is very doubtful in the colour- 
suppressed case, it should work much better for the colour- allowed channels because 
of colour transparency JT^]. The latter have furthermore rather tight restrictions from 
the Heavy Quark Effective Theory (HQET) |25|] for the form factors describing the 



"factorized" hadronic matrix elements of the relevant four-quark operators [24], [26 

Our paper is organized as follows: in Section [D] we calculate the transition matrix 
elements and observables of the angular distributions by using an appropriate low-energy 
effective Hamiltonian. There we also give estimates for these observables, allowing a com- 
parison with experimental data. The efficient experimental determination of these ob- 



servables is the subject of Section |TJ, where we shall discuss the angular moment analysis. 
Sections |TV] and [V] are devoted to the angular correlations in the colour-suppressed decays 
B s -> J/ijxj), B -> J/ipK* and the colour-allowed decays B s -> D*+D*-, B -> D* + D*, 
respectively. There we give the time evolutions of the angular distributions, appropriate 
weighting functions, and discuss CP-violating effects. Finally in Section [V| the main 
results are summarized. 



II Transition matrix elements and observables 

Before we present an efficient method for extracting the observables of the angular dis- 
tributions from experimental data - the angular moment analysis - let us discuss in this 
section how these observables are calculated and what orders of magnitude we expect 
for them. 

II-A General aspects 

In order to calculate the decay amplitudes of the b — > sec transitions considered in this 
paper, we use an appropriate low-energy effective Hamiltonian, which has the following 
structure: 



( 10 

E M s) $Ci(aO + <&cm + E QkCM 

j=u,c K k=3 



(2) 



(s) 

Here the quantities A] = Vj S V* b denote CKM factors, 



Ql — (c a S /3 )v-A (fr/3C a )v-A > Q2 — ( c aS a )v-A (ft^C^jy-A 

Q1 = (u a S{3)v-A (&/3M a )v~A , Q\ = (M Q S a )v_ A (bpU^v-A (3) 

are "current-current" operators, 

Q3 = (b a S a )y- A E (Wp)v-A , Qa = (baSp)v-A E (Wa)v-A 
q=u,d,s,c,b q=u,d,s,c,b 

Q5 = (b a s a )v-A E (QpQph+A , Qe = (&aS /3 ) V -A E (Wah+A (4) 

q=u,d,s,c,b q=u,d,s,c,b 
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describe QCD penguins, while the operators 



3 3 - 

q=u,d,s,c,b q=u,d,s,c,b 

3 3 - 

Q9 = 2 ( b aS a )v-A e 9 (Wph-A , QlO = 2 ( & a S /3)v-A 6 Q (QpQa)v-A 

" q=u,d,s,c,b ~ q=u,d,s,c,b 

( 5 ) 

are "electroweak" penguin operators. Here V ± A corresponds to 7^(1 ± 75) quark cur- 
rents, Greek indices are associated with the S77(3)c quark-colour, and the quantities 
e q arising in the expressions for the electroweak penguin operators label the electrical 
quark charges. Nowadays, the Wilson coefficient functions Ck(fi) of the low-energy ef- 
fective Hamiltonian 7i e s, where fi = 0(mb) denotes the usual renormalization scale, are 



known beyond the leading logarithmic approximation [27 



Since A^ is suppressed with respect to A^ by a CKM factor \ 2 Rb, where A = 0.22 



is the Wolfenstein parameter [10] and 



1 


v ub \ 


A 


\Vcb\ 



is constrained by present experimental data to lie within the range Rf, = 0.36 ± 0.08 
|28| . |29| , and since furthermore the current-current operators Q", Q2 may contribute only 
through penguin-like matrix elements to b — > sec modes, the corresponding transition 
amplitudes are dominated to an excellent approximation by the contribution proportional 
to A^ (for a detailed discussion, see |30[| ). In the penguin operators, we neglect the 
parts of flavour structure different from (cc) ihs) . Then the number of relevant operators 
reduces from ten to four and the structure of the decay amplitude simplifies considerably. 

In order to implement the factorization hypothesis by factorizing the hadronic matrix 
elements of the four-quark operators Qk into hadronic matrix elements of quark currents, 
we have to perform suitable Fierz transformations of the operator basis specified in (||)- 
(||). Beyond the leading logarithmic approximation one has to be very careful in perform- 
ing such Fierz transformations, as the Wilson coefficients depend both on the form of the 



chosen operator basis and on the applied renormalization scheme ||27|| . Since we do not 
use any specific Wilson coefficients to obtain numerical estimates in this paper, we may 
perform such Fierz transformations and will use a tilde (~) to indicate Fierz-transformed 
operators. For a discussion of the renormalization-scheme dependences arising beyond 
the leading logarithmic approximation and their consistent cancellation in the physi- 
cal transition amplitudes through certain one-loop matrix elements at fi = 0(mb), the 
reader is referred to Ref. |TT . 



Let us, in the following two subsections, investigate the structure of the hadronic 
matrix elements of the low-energy effective Hamiltonian [Eq. (0)] for the exclusive colour- 
suppressed and colour- allowed decays B s — > J/ip (ft, B — >• J/ipK* and B s — > D^ + D*~, 
B — > D* + D , respectively. 
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II-B Colour-suppressed decays 

If we perform a Fierz transformation of the current-current operators specified in (^) , the 
decay amplitude for B q — > J/ipV ((q, V) G {(s, 0); (d, K*°); (u, K* + )}) can be written in 
the following form: 

(j/4j(\)v(x)\n cS \B q ) = ^v cs v; b (7) 



x 



Cf( / .)(J/^(A)y(A)|QJ( / ,)|5 g )+Cf oct ( / z)(J/^(A)y(A)|Qtoct(/^)l^) 



+Cf( / i)(J/^(A)V(A)|g^(^)|^)+Cf oct ( / x)(J/^(A)V(A)|g c 5>oct ( /U )|^) 

where A denotes the helicities of the final-state vector mesons and the "effective" Wilson 
coefficient functions are given by 

cfOi) = cm + \cm + cm + 1 -cm + CM + ic 10 (//) (8) 

ct%M = 2 [cm + CM + CM] (9) 

cf(fi) =cm + \cM + cm + \cM (io) 

ci%M = 2[cm + cm\- (ii) 

The /^-dependence of these Wilson coefficients is cancelled by that of the hadronic matrix 
elements appearing in Eq. (|7|). In deriving the transition matrix element in Eq. (^), we 
have used the relations 

Q c 2 = 1q$ + 2Ql oct (12) 
Ql = \Ql + 2QU (13) 



with 



and 



Qi = (c*c*)v-a (h s p) v _ A (I 4 ) 

Qloct = faZp<*)y_ A (M>)v-A ( 15 ) 

Ql = (c Q c a ) v+A (bpSp) y _ A (16) 

<?U = iS^cp) (M>) V _ A • (17) 



' V+A 



Here the 3x3 matrices T a are the SU(3)c generators, normalized to tr(T a T ) = 5/2. 
As we will see below, the form of the Fierz-transformed operators given above is better 
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suited to analyse the B q — > J/ipV decays since the J /if) is related to the (cc) pieces. 
The penguin contributions to Cf s (fj,) and Cf s oct (fi) are at most 0(10%) and 0(1%), 
respectively, as can be estimated from the values of their Wilson coefficients |27j . 

If one assumes that J/ip emerges from the vector parts of the (cc)v±a quark currents 
appearing in the operators in Eqs. (|T^ ) - ([T7|) , the matrix elements of Qi( jOC t) and Q5( jOC t) 
will be equal and the decay amplitude Eq. (|7|) can be simplified considerably. Moreover, 
within the framework of naive factorization, we obtain (analogous for Q\ and Q5 OCt ): 



JMA)V(A) Q\ 
JMA)V(A) Qi oct 



B, 



B, 



m(\)\(cc) y _ A \0)(v(X)\(bs) 
(J/^(X)\(cT a c) v _ A \0) (V(X) 



V-A 

bT a s 



B. 



V-A 



B. 



(18) 
(19) 



where summation over colour- indices is understood implicitly. Consequently, since J/ip 
is a colour- singlet state, the factorized hadronic matrix elements of the colour-octet 
operators given in Eq. ( |i~9"|) vanish. 



II-C Colour-allowed decays 

In the case of the colour-allowed decays B q — > D* + D* q (q 6 {u,d,s}), the transition 
amplitude can be written in a way that is completely analogous to Eq. (|7|): 



(D: + (X)Dl(X)\H cS \B q ) = ^V cs V: b 



(20) 



x 



Cf^){D: + {X)D^X)\Q^)\B q )+Cll^ 

+Cf(ri(D: + (X)^(X)\Q^)\B q )+C^M ■ 
The corresponding effective Wilson coefficient functions are, however, very different: 



cf(fi) = -cm + cm + -cm + cm + ~cm + cm 



cf M = ncM + cM + cM\ 



cfM = ~cm + CM + ~c 7 (») + cm 



i 



neff 



3 v ' vr/ 3 

' 6l octM =2[cm + cM\- 

In deriving Eq. (^), we have used the relations 



Ql ~ + 2<52,oct 



Ql 



1 X 



Ql + 2Q C 



6, oct 



(21) 
(22) 
(23) 
(24) 

(25) 
(26) 
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with 



Q c 2 = (c a s a ) v _ A (bpCp) y A 



V-A 



and 



Qe,oct 



-2 (c a Ls a ) (bpRcp) 

-2 (c a LTZp S p) (b 7 RTy s 



(27) 
(28) 

(29) 
(30) 



Here L and R correspond to the Dirac structures 1 — 75 and 1 + 75, respectively. The 
D* + meson emerges from the (cs) pieces of these operators. Since it is a vector meson, 
we have 

[D* s + \c a Ls a \0) = 0, (31) 



and hence the factorized matrix element of Q\ vanishes. As in Sec. il-B , the hadronic ma- 
trix elements of the colour-octet operators vanish within the factorization approximation 
because of their colour-structure. 



II-D Observables of the angular distributions 

The hadronic matrix element of a generic four-quark operator Q between the state vectors 
(Vi(A)V2(A)| and \B q ) has the following general Lorentz-decomposition ]5], ||: 



<yi{\)V 2 {\)\Q\B q ) 
ev llAt (A)*ey 2il ,(A)* 



U V 1 



m Vl m V2 



-e^p Vl ,aPv 2 ,p 



(32) 



where the symbols e(A) denote the polarization vectors of the final-state vector mesons 
V\ and V2. A similar parametrization can be employed to express the transition matrix 
elements [Eqs. (0) and fl20Dl, yielding 



a = ^V cs V cb 



Cf^)A\ + Cf + ^)A\ +4 (33) 
+ Cf{p)Af{n) +(^{p)KM +C%MAf +4 (/i) +C^ 4jOCt (/i)A n + 4 ,oct(^)" 



Cf^)B\ + C^)B\ +4 



cff 



b = T2 VM 



(34) 
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^-V V* 

V2 VcsVcb 



Cf{»)Cl + C^)Cl +A 



efi 



(35) 



+ cf {fief (ai) + c5 ct (^)cfoct(^) +c?U») c £M +cf + *,oMc? +4 , M 



inf/ 



7eflf 



mf 



cff 



tnf 



->cff 



mf 



where the index i distinguishes between colour-suppressed {i = 1) and colour- allowed 
(i = 2) decays and "f" and "nf" correspond to "factorized" and "non-factorized" matrix 
elements, respectively. Note that the factorized amplitudes do not depend on the renor- 
malization scale /i. Since the Wilson coefficients depend on this scale, this already signals 
the need for non-factorizable contributions to cancel the /i-dependence in Eqs. 



V1V2 in terms of linear 



(see e.g. Ref. for a further discussion of that point). 

In the following sections we will analyse the decays B q 
polarization states. The corresponding decay amplitudes take the form p3| [33 

A {t) 



A(B q (t) -> V 1 V 2 ) 



x 



*L *L 

e v 2 



A l{ (t)e 



Vi e v 2 



2-iA ± (t)e* Vl xe* V2 -p v jV2, (36) 



where x = pv 1 • Pv 2 / i m v 1 m v 2 ) an d Pv 2 is the unit vector along the direction of motion of 
V2 in the rest frame of V\. Here the time dependences originate from B q -B q mixing. In 
our notation, an unmixed B q meson is present at t = 0. 

The linear polarization amplitudes at t 
terms of a, b and c as follows |03| : 



A o (0) 
A\\(0) 

Ajl(0) 



— xa 
V2a 



defined by Eq. (|36|) can be expressed in 
- (x 2 - l)b 



(37) 



2(x 2 



1) c. 



At time t 
IA)(0)|, \A 



0, the angular distributions for B q — > V1V2 depend on the observables 



and on the two phases 5\ = Arg An {0)*A±(0) 



and 5o 



Arg [y4o(0)*v4j_(0)], which are CP-conserving strong phases that are (mod ir) in the 
absence of final-state interactions (probably not a justifiable assumption for the colour- 
suppressed modes). Quantitative estimates for these observables will be given in the 
following subsection. 



II-E Factorization tests and estimates of observables 

While the non-factorizable contributions to a, b and c cannot be calculated at present, the 
evaluation of the factorizable contributions is straightforward. Without yet going into the 
details of which form factors to employ, the naive factorization assumption yields many 
testable consequences. For example, time-reversal invariance forces the form factors 
parametrizing quark currents to be all relatively real. Consequently, naive factorization 
predicts the same strong phase (mod 7r) for the three amplitudes A (0), A\\(0), A±(0). It 
therefore predicts vanishing values of the two observables |L7|, [L9|, |20| 

Im L4*(0)A ± (0)] = (38) 
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Im [Al(0)A ± (0)} = 0, (39) 

and the equality 

Re[^(0)A|,(0)] = ±|A (0)A||(0)|. (40) 

The breakdown of the naive factorization assumption is unequivocally proved if any of 
the three equations (|38l)-(f40|) is not satisfied. Detailed comparisons of polarization am- 
plitudes in non-leptonic and semi-leptonic decays test additional implications of the naive 
factorization assumption. The phenomenology of detailed studies of the full non-trivial 
angular distributions is thus much richer than the single factorization test available for a 
pseudoscalar decaying into two pseudoscalars |18[. While the above equations represent 



general tests of the factorization assumption, it is also useful to examine the predictions 
for the observables of the angular distributions for various form factor ansatze. 

II-E.l The colour-suppressed decays B q — > J/ipV 

The factorized amplitudes for B q — > J/ipV with (q,V) G {(s,<f>); (d,K*°); (u,K* + )} are 
given by fl6], ^T], ^2] : 



A[ = - fj/fmj/ f (m Bq +m v )Af qV (m 2 j^) = A[ 

B 1 2 — - — A 2 (m^) - B 5 (41) 

C{ = 2 fj/ ^ mV V^(m^) = Cl 
1 m Bq +m v v J/i,J 5 ' 

where we have used the notation of Bauer, Stech and Wirbel for the form factors Af qV (q 2 ) 
and V BqV (q 2 ) of quark currents |TB[. The parameter fj/^p denotes the J/ip decay constant, 
which can be determined from the J/ip e + e~ rate, yielding fj/^ = 395 MeV. 

At present, several methods for obtaining the form factors A 1 (mj / ,^), A 2 (mj/^) and 
K(mm) for the B — > K* case are on the market. Using SU(3) flavour symmetry of 
strong interactions, the B — > K* form factors can be related to the B s — > <fi case. In 
Table |] we have collected the form factors proposed by several authors [[16], |34|, ^ 



and have moreover given the corresponding predictions for the ratios of observables of 
the angular distributions. These ratios should suffer less from unknown SU (3)-breaking 
corrections than the observables themselves. Note that these ratios are independent of 
the Wilson coefficients within the factorization approach. 
The quantity 

t (r\\ i a lc\\\2 

(42) 



r (o) 


A)(0)| 2 


r (o) + r T (o) 


A,(o)| 2 + 


^||(0)| 2 + 


^4±(0)| 2 



describes the ratio of the longitudinal to the total rate at t = 0. Although CDF |36 



claims to have measured this quantity, from their untagged data sample, to be 0.56 ± 
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Observable 


BSW 16 


Soares |34j| 


Cheng H 




0.46 


0.42 


0.41 




0.46 


0.43 


0.36 


V»«'{mfo) 


0.55 


1.08 


0.72 


|A||(0)|/|A o (0)| 


0.81 (0.77) 


0.82 (0.78) 


0.75 (0.70) 


\A ± (0)\/\A (0)\ 


0.41 (0.40) 


0.89 (0.88) 


0.55 (0.54) 


r (o)/(r (o) + r r (o)) 


0.55 (0.57) 


0.40 (0.42) 


0.54 (0.56) 


5i 


71 


71 


71 


5 2 












Table 1: Predictions for form factors and B s — > J/ipcj) (B — > J/ipK*) observables. 

0.21(stat.)±o;[)4(syst.), their claim is valid only if the CP-odd component of B s — > J/ip<t> 
is negligible, or if the lifetime difference Ar can be ignored. 

The 2 nd -4 th columns of Table [l] are calculated within the framework of naive fac- 
torization, i.e. we have inserted Eq. (|4"ID into Eqs. (|33|)-(|3"5|) and have omitted the "nf" 
terms in order to calculate the amplitudes in Eq. (|37|). 

The form factors given by Soares [Bl] are obtained from D — > K^Hvi data by using 



heavy-quark symmety relations [3^| and assuming the monopole momentum-transfer 
dependence of the BSW model [16]. Some more form-factor models and their predictions 
are discussed in . Note that the small difference between the B s — > J/ip <ft and 
B — > J/ipK* results in Table [I] is related to phase-space effects and not to any SU(3)- 
breaking effects in the corresponding hadronic matrix elements. 

Looking at Table [l], we observe that the "factorized" predictions for |A|i(0)|/|Ao(0)| 
are rather stable (~ 0.8), while |Aj_(0)|/|A)(0)| depends strongly on the method used 
for obtaining the form factors. A common feature of all results is 5\ = 71 and 62 = 0. 
Therefore a measurement of non-trivial phases 5i and 62 would imply the presence of 
strong final-state interactions and non-factorizable contributions. 

Whereas the use of the factorization assumption is very questionable in the case of 
the channels B s — ► J/ip <fr and B — > J/ipK*, flavour SU(3) symmetry is probably a good 
working assumption. Thus all the hadronization dynamics of the B s — > J/ip4> decay, 
such as the phases 61 and 5 2 and magnitudes of the amplitudes 

A)(0), A|,(0), A ± (0), (43) 

can be obtained from the B — > J /ipK* modes.Q This approach may be helpful to extract 
the CKM phase 0ckm (see Eq. ([]])), as we will see below. 

1 Although those SU(3) relations are mostly trivial, one subtlety due to quantum-coherence must be 
emphasized. Because of the SU (3) relations in the unmixed amplitudes 

Af(B s — > J/ip(f>) = Af(B — > J/ipK*) , where / = 0,||, ± , 

the magnitudes of the amplitudes for or Bf decays into CP-even or CP-odd J /ip(j) final-state 
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The factorization assumption should work much better for the transitions B s — > 
D* + D*~ and B — > D* + D . Therefore the results presented in the following subsection 
should be more reliable than those summarized in Table |I|. 



II-E.2 The colour-allowed decays B q -> D* + D* q 

Using again the same notation as Ref. [fL6|, we get the following "factorized" results for 
the modes B q -> D* s + T3* (q G {«, d, s}) H, g| |f 



^2 = -fD* s m D *(m Bq +m D *)A 1 



Bn D 



m 



C 2 f = 2 



f D *ml t m D * aBqD 
-i 

fD*rn 2 Dt m D: 



"^b + m D * 



A = 



B'=0 



a = o. 



(44) 



The parameter fp* is the D* decay constant. The spin symmetry of HQET implies fr>* ~ 
fo s - A recent compilation of measurements of fn s from D s — > /iI7 gives (241 ± 21 ± 30) 



MeV 39 



In the case of B q — > D q transitions we have rather tight restrictions from HQET (for 
reviews, see for example ||25| ) for the corresponding form factors. The following ratios 
turn out to be useful to implement these HQET constraints [fl0|j : 



R2H 



1 



m Bq + m D , 



m Bq + m D * 



Af D *{q 2 ) 

A$ D *(q 2 ) 
Af D *(q 2 Y 



where Ri(w) and R2{w) are defined in such a way that we have 

RAw) = R 2 (w) = 1 



(45) 



(46) 



(47) 



for all values of w in the strict heavy-quark limit. The kinematical variable w is defined 
by 



w 



m B q + m D q ~ (f 

2m B m D « 



(48) 



configurations, respectively, are a factor of \[2 larger than their corresponding B — * J/tpK* ones. [Here 
the K* is seen in a flavour-specific mode. If K* is neutral and is observed as ir°Ks, quantum coherence 
in B° — B° must also be taken into account.] If the CP-even processes dominate, then 

r(i? s L -» j/u) « 2r(s -> j/i>k*) . 

Studies of B s versus B production fractions can thus be undertaken, since the lifetimes will be precisely 
known. 
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Observable 


BSW 


HQET c +„:„ f 

v strict; 


^ sym. .- DrGaK. 


A» u *{m 2 D ,) 


0.72 


0.70 (0.68) 


0.70 (0.68) 




0.76 


0.76 (0.75) 


0.54 (0.53) 




0.79 


0.76 (0.75) 


0.89 (0.88) 


|An(0)|/|A Q (0)| 


0.90 (0.90) 


0.91 (0.91) 


0.81 (0.80) 


\A ± (0)\/\A (0)\ 


0.32 (0.32) 


0.32 (0.33) 


0.33 (0.34) 


r (o)/(r (o)+r T (o)) 


0.52 (0.52) 


0.52 (0.52) 


0.57 (0.57) 


Si 


7T 


7T 


7T 


s 2 












Table 2: Predictions for form factors and B s — > D* + D* (B — > D* + D*) observables. 



The value of the momentum transfer q 2 relevant for Eq. 
factor Af D * (q 2 ) is usually written as 



2y/mB q m D * 



m Bq + m D * 



is q 



h Al (w) 



m 



Df 



The form 



(49) 



where h Al (w) corresponds to the Isgur-Wise function in the strict heavy-quark limit and 
can be written as 



h Al (w) = 1 - p 2 Ai (w - 1) + 0((w - 1) 



(50) 



The current status of the normalization ^(1) and of the "slope parameter" p\ has been 
summarized recently by Neubert in Ref. |4(|. The form factor Af D *(q 2 ) is protected 
by Luke's theorem E1J against I/tuq corrections at zero recoil. The other form factors 
A-2 D * {(f) and V BD * (q 2 ) are not protected by this theorem. From calculations based on 
HQET one expects a rather weak dependence of R\(w) and R2{w) on w and therefore 
uses 



R^w) = R^l + O^w- 1)] 
R 2 {w) = R 2 [l + 0(w- 1)]. 



(51) 
(52) 



In our analysis we will neglect the u>-dependence completely. 

Following these lines we have calculated the results for the form factors and ratios 
of observables, which should receive smaller S'?7(3)-breaking corrections than the ob- 
servables themselves, summarized in Table 0. For completeness we have also given the 
results obtained by applying the BSW model []IB| in the 2 nd column. In order to calculate 
the 3 rd and 4 th columns, we have used .F(l) = 0.91 and p\ = 0.91 
denoted by HQET strict and HQET sym _ break correspond to R± = R 2 - 
R 2 = 0.71, respectively, where we have employed the results by Neubert p(| to take into 
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The columns 
1 and R 1 = 1.18, 
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account HQET symmetry-breaking corrections. Within the factorization approximation 
we obtain the following simple expressions for A\\(0)/Aq(Q) and A±(0) / A (0) in terms of 
the HQET parameters: 




-i -i 



A)(0) 



i?2(w) — x 



(53) 



A ± (0) 



4,(0) 



) 



Ri(w) 



A (oy 



(54) 



where the kinematical variable 



x = 




(55) 



2m D *m D * 



has been defined after Eq. (|36|) . 

If we compare Table § with Table [l], we note that the results for the observables 
depend much less on the way of obtaining the form factors. Also the "old" BSW model 
is in rather good agreement with the HQET predictions, which is quite remarkable. 
Therefore the results given in Table | are more reliable than those collected in Table |l|. In 
this respect it is also important to note that the non-factorizable contributions appearing 
in Eqs. (|3^)-(|35|) should play a minor role for the colour- allowed decay class and that 
5\ = 7r and <5 2 = is expected to hold on rather solid ground. Because of the latter 
feature, CP-violating effects arising in untagged B s — > D* + D*~ data samples should be 
a promising way to extract the weak phase 0ckm (see Eq. ([[])), as has been outlined in 
detail in Ref. §. 

Ill The angular- moment analysis 

The main focus of this section is the efficient determination of the observables discussed 
in Sect. [I I^D] and [11- E| . This can be accomplished by an angular-moment analysis ffl]. 
In this approach, the observed data are weighted by judiciously chosen functions, which 
project out any desired observable. Whereas Ref. |1| determines the moments for a 
few choice angular distributions, using spherical harmonics, this paper indicates how to 
determine suitable weighting functions for all kinds of angular distributions, using only 
orthogonality arguments (without invoking spherical harmonics). 

Let us denote the angular distribution of a given decay by 



where a represents all the parameters that are independent of the kinematics, which is 
described by certain decay angles. In general, the physical process involves an arbitrary 



f(e,a;t) =$>«(«;*)<?« (9), 



(56) 
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number of such angles denoted generically by 0. For the examples considered in this 
article, % runs from 1 to 6, and we have 



a = {T Hl T L , Am, L4 (0)|, |A|,(0)|, |4l(0)|, 8 U 5 2 , 0ckm}- 

All the quantities of interest are encoded in the time evolution of the observables (a; t). 
In the following discussion the a- and t-dependence of the &W's is implicit wherever not 
explicitly stated. 

The usual method for extracting ftW's is to use an unbinned maximum likelihood fit 
Q . Performing such a fit for a given quantity requires some idea of the values of the other 
quantities. When one deals with limited statistics, one may want to exploit alternative 
methods that completely decouple the extraction of one observable from all the others. 
Luckily such a method exists, the angular-moment analysis. 

If we can find a weighting function u;W(G) for each i such that 

/[^e]«,W(e)^)(e) = < y tt = |j^jf * (57) 

then the feW's can be obtained directly from 

&« = J[De]w®(e)f expt (e) . (58) 

Here [DQ] denotes the appropriate measure for integrating over all angles 9, and / ex pt(@) 
denotes the observed full angular distribution. For a small number of events (N), the 
form of the function / ex pt(@) will not be known, but only the values for will be known 
for each event. In that case the above equation reduces to 

&W = i £ u;W(6) . (59) 

events 

These &W's can then be used directly for studying their (a;t) dependence. 

That i«w's can always be found for any angular distribution follows from the linear 
independence of the gfW's (they have to be independent for the angular distribution to 
be legitimate). The vector space V& spanned by all g^'s for j k is a proper subspace 
of the vector space V spanned by all g^'s. Then there exists a one-dimensional vector 
space W such that V = Vfc © W and Vfc _L W. Here the scalar product is defined as 
Vi ■ v 2 = J[DQ]v 1(6)^2(0); is then the element of W with proper magnitude. 

For a given set of g'^'s, the choice of tyW's need not be unique. We can always 
take any vector space V' D V and the corresponding projection space W" such that 
V = V k ® W and V fc JL W. Then any w (fc) G W with w (k) ■ g ik) = 1 will serve our 
purpose. 

We now indicate an explicit procedure for finding a set of weighting functions ap- 
plicable to any given angular distribution. For a theoretical angular distribution of the 
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form f = b®g® (the dependence on angles and time is implicit), 

i=i 



^ W = EW J) (60) 

i=i 

is a proper weighting function, where the n 2 unknowns A^- are solutions of the n 2 simul- 
taneous equations 

Sik = x v J [DO}g U) g {k) . (61) 

The existence of such a solution follows from the vector-space arguments given earlier. 
The iuw's need not be restricted to the vector space spanned by the n vectors g&% in 
which case the unknowns Ay will be underdetermined and more than one set of w^'s 
will serve our purpose. 

It is crucial to observe that the weighting functions w^ 1 ' depend only on the angular 
terms and not on the values of the observables The implication is that no matter how 
complicated the detailed angular distribution, there always exists an angular weighting, 
which projects out the desired observables alone. We therefore recommend the use of 
moments whenever one wishes to extract observables from measured angular distribu- 
tions, such as in weak decays of baryons or pseudoscalars [P — > V£u, Xjiv, VV, etc.], 
or strong and electromagnetic decays |J . The utility of this approach cannot be overem- 
phasized. For instance, the moment analysis allows the study of the q 2 dependence of 
each of the observables separately in the process P — > Viv. This could prove useful for 
the extraction of form factors and the determination of CKM elements, e.g. V c b and V^&. 

We note that there exist many legitimate choices of weighting functions. The opti- 
mal choice depends on the numerical values of the observables [IJ and on the detector 
configuration. 

IV The angular distribution of the colour-suppressed 

decays B s — > J/ip (f> and B — > J/ipK* 

In this section we give the angular distribution of the decays B s — > J/ip <fi and B — > 
J/ijjK*, their time-dependences and appropriate weighting functions. 

IV-A The decay B s -> J/^(-> K+K ) 

An analysis of this process has been performed in |33] in terms of linear polarization 
states of the final-state vector mesons. The corresponding decay amplitude has the same 
form as Eq. fl5E|). Since the amplitudes Aqh and A± are related to CP-even and CP-odd 
final-state configurations, respectively, they differ in time evolution as well as angular 
distribution. The angular distribution can be used to separate these components and 
their time evolution can be studied individually. 
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The differential decay rate at time t as a function of a generic variable x will be 
denoted by 

dT(t,x) ^ 1 d 2 N(t) 

dx N{t) dxdt ' 1 ' 

Consequently the normalized number of decays in the intervals [t, t + At] and [x, x + Ax] 
is given by 

dT(t,x) A 1 d 2 N(t) x A , N 

V Ax At = -——^AxAt. 63 

IV-B Tagged decays 

In the case of 5 S — > J/^ 0, the three-angle distribution for the decay of an initially 



present (i.e. tagged) _B S meson takes the form |33 



d 3 T[B s (t) -> JM(-> Z + rW-> K+K-)} 9 



2|A (t)| 2 cos 2 ^(l - sin 2 ^cosV) 



dcosO dip dcosip 32tt 
+ sin 2 ^{|A||(t)| 2 (l -sin 2 #sinV) + |A ± (t)| 2 sin 2 9 - Im (Aj|(t)A ± (t)) sin 29 sin y?} 

+ ^sin2^{ Re (A*(t)A||(t))sin 2 ^sin2^+ Im (A* (t)A x (t)) sin2#cos^ } 1 . (64) 
a/2 J 



Throughout this section we will apply the same conventions as in Ref. |33[], i.e. <p moves 
in x direction in the J/ip rest frame, the z axis is perpendicular to the decay plane of 
4> — ► K + K~ , and p y (ir + ) > 0. The coordinates (6*, <p) describe the decay direction of l + 
in the J/ip rest frame and ip is the angle made by p(K + ) with the x axis in the (j) rest 
frame. With this convention, 



--Pd>(p^-p A -+) 
--P4>(P0-Pjr+)l 

sin 6* cos y? = • x, sin 9 sin (p = • y , cos 9 = p£+ ■ z . (65) 



y = I ^ ^-y-r z = x X y 
J Ip k +-P0(p^'Pk+)I 



Here, the bold-face characters represent unit 3- vectors and everything is measured in the 
rest frame of J/if). Also 

cos^ = -p' K + ■ p' J/v „ (66) 

where the primed quantities are unit vectors measured in the rest frame of 0. 

The time dependence of the right-hand side of Eq. (64) can be read off from Table [| 
where Am = m# — mi > is the mass difference of the mass eigenstates Bf and 
of the B s system and T = (r# + Tl)/2 denotes their average decay width. The 
phases 5\ = Axg(A\\ (0)*A±(0)) and 5 2 = Arg(A (0)* A±(0)) are CP- conserving strong 
phases. In the absence of final-state interactions - probably not a justifiable assumption 
for B s — > J j ipcf) - they are expected to be (mod tt) . 

On the other hand, the quantity 5(f> = 0ckm (see Eq. ([!])) is a CP-violating weak 
phase, which is introduced through interference effects between B s -B s mixing and decay 
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Observable 


Time evolution 


\Mt)\ 2 
I^iiWI 2 
\Mt)\* 


l^o(0)| 2 
l^||(0)| 2 
l^(0)| 2 


e Tht — e rt sin(Amt)5<j) 
e~ Vht — e~ rt sin(Amt)<50 
e -r H t + e -Fi s i n (Amt)^ 




Re(A*(t)A {l (t)) 
lm(Al(t)A ± (t)) 
lm(A*(t)A ± (t)) 


|A (0)||A||(0)| cos(<5 2 - 5 X ) [e- r ^ - e~ Tt sin( Amt) 5 cj)\ 
|A||(0)||A ± (0)| [e-^sin^! - Amt) + \ (e~ rBt - e~ v A cos(<5i)50] 
|A)(0)||A ± (0)| e- r *sin(5 2 - Amt) + \ [e^ - e~ r ^) cos(<5 2 )<50 



Table 3: Time evolution of the decay B s — > J/if>{— > Z + / )</>(—► K + K ) of an initially 
(i.e. at t = 0) pure 5 S meson. 



Observable 


Time evolution 


W*)l a 

l4L(t)l 2 




4>(0)| a 
^||(0)| 2 


e~ Tht + e~ r * sin(Amt)<50 
e~ Tht + e~ rt sm(Amt)5(fi 
e -r H t _ g-rt s in(Amt)5(/) 




Re(A* (t)A {l (t)) 
lm(Al(t)A ± (t)) 
lm(A;(t)A ± (t)) 




A (0)P||(0)| cos(<5 2 - Si) [e~ TLt + e~ Tt sm(Amt)5(j)\ 
-|A||(0)||Aj.(0)| [e" r * sin(5i - Amt) - \ (e~ FHt - e~ T ^) cos(5i)<tyl 
-\A {0)\\A ± {0)\ 'e- Tt sm(5 2 - Amt) - \ (e"^' - e~ r ^j cos{5 2 )5<j> 



Table 4: Time evolution of the decay B s — > J/if>(— > Z + Z )4 > {~^ K + K ) of an initially 
(i.e. at t = 0) pure £? s meson. 



processes. It can be expressed in terms of elements of the CKM matrix || [12| as 

^ VcsVcb 



exp(i5(f>) 



v t *v tb v cs v; b 



(67) 



and is very small, as can be seen easily by applying the Wolfenstein expansion of the 
CKM matrix flltl 



At leading order in this expansion 



vanishes. However, taking into 
account higher-order terms (for a treatment of such terms, see e.g. Ref. |4^| ) gives a 



non- vanishing result |13 



2X 2 V = 0(0.03). 



sin Or 



(68) 



0.22 is 



Consequently 8(f) measures simply the CKM parameter 77. Note that A 
related to the Cabibbo angle. Useful expressions for 8(f) can be found in Ref. fL3f , where 
the following relation has been derived: 



2X 2 R b sin 7. 



(69) 
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Here 7 is one angle of the "usual" unitarity triangle [44|. Consequently, if the CKM- 
parameter Rf, (defined by Eq. (|])) is used as an input, 5(f) allows a determination of 7. 
That input allows, however, also the determination of rj (or 7) from the mixing-induced 
CP asymmetry of — > J/ipKs measuring sin 2/3, where j3 denotes another angle of 
the unitarity triangle f44|| . If one compares these two results for 77 (or 7) obtained from 
B s and B^ modes, respectively, a test of whether the B s -B s and B^-B^ mixing phases 
are described by the Standard Model, or receive additional contributions from physics 
beyond that model can be performed. Needless to note, a measurement of a value of 
5(f) much larger than the Standard Model expectation of (9(0.03) would anyway be a 
striking signal for new physics in B s -B s mixing. 

An interesting interpretation of 5(f> has been given in Ref. jRJ. There it was shown 
that 5(f) is related to one angle in a rather squashed (and therefore "unpopular" ) unitarity 
triangle. Note that terms of 0(5(f> 2 ) have been neglected in Table [| 

The angular distribution for an initially present B s meson is given by 



d 3 T[B s (t) -> J/V>(-> K + K-)} 9 



dcosO d(p d cos if) 



oc 



32tt 



2|A (t)| 2 cos 2 ^(l -sin 2 flcosV) 



+ sin 2 ^{|A||(t)| 2 (l - sin 2 #sin 2 ^) + \A ± (t)\ 2 sin 2 9 - Im (A\(t)A ± (t)) sin2#sin<^} 



+ -=sin2^{ Re (A;(t)A||(t))sin 2 #sin2^+ Im (A* (t)A ± (t)) sin 28 cos p } 
a/2 



(70) 



where the angles are again defined by Eqs. ( p5|) and ( pq) . The time dependence of 
this rate can be obtained easily with the help of Table [|, where terms of 0(5(p 2 ) have 
been neglected, as in Table |^. In calculating Tables [| and |] we have used the fact 
that B s — > J/ip(f> (and B s — > J/ifxp) is dominated to excellent accuracy by a single 
weak amplitude, as we have seen in Section |I[ Therefore we have to deal only with 
mixing-induced CP violation and there is no direct CP violation, i.e. |Ao(0)| = |Ao(0)|, 
|A||(0)| = |Aj|(0)|and|Z ± (0)| = |A ± (0)|. 

It is important to note that the mass difference Am can be extracted from time- 
dependent analyses of tagged B s — ► J/ifxf) data samples Previous experimental 
feasibility investigations for the extraction of Am focused entirely on tagged flavour- 
specific modes of B s mesons [pEE|, E7|. 



IV-C Untagged decays 

Combining Tables |3] and [|, we find that the time evolution of the untagged data sample 
for / = J/V>(-> K+K~) is given by 



d 3 T[f(t)} ^ 9 



dcosO dp dcosip 16w 



2|A (0)|V ri ' cos 2 ip(l - sin 2 cos 2 p) 



+ sin" ^{| Ay (0) \ 2 e~ TLt (1 - sin" sin" y>) + | A±(0) \ 2 e' L Ht sin" 0} 
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+-J= sin 2ip {| A o (0) \\A\\ (0) | cos(5 2 - 5i)e~ rL * sin 2 9 sin 2y?} 
+ 1 -j= | A (0) 1 1 A± (0) | cos 5 2 sin 2^ sin 29 cos <p 

- |A||(0)p ± (0)|cos5isin 2 ^sin2^sin^|^ (e- rff '-e~ rL *)(50 . (71) 

Remarkably the time dependence of the untagged rate does not depend on the mass 
difference Am. This feature has been discussed within a more general framework in 



Ref. |jI4}| . Consequently, whereas I\ and T# can be determined from the untagged 
data sample, the extraction of (Am)g s requires tagging. As has already been pointed 
out in ||, because of the lifetime difference (AT)b s , the untagged decay rate [Eq. flTID l 
develops an interesting contribution for t > 0, which is proportional to the CP-violating 
weak phase 5<f). It originates from the imaginary parts of the interference terms between 
A±(t) (A±(t)) and Afl(t) (Aj(t)), A* (t) (A* (t)). li AT = T H - T L is in fact sizeable, we 
are optimistic that it will be possible to measure this effect. 

IV-D A closer look at the one-angle distribution 

The full three-angle distributions for tagged and untagged B s — > J/if>(—> Z + Z~)0(— > 
K + K~) decays discussed in the previous subsections are quite complicated. A much 
simpler case arises if we integrate out the two decay angles (p and if) in ([64]), leading to 
the following one- angle distribution |33 : 



oc (\A (t)\ 2 + L4||(t)| 2 )f(l + cos 2 9) + \A ± (t)\ 2 hinH. (72) 
a cos 9 8 4 

Let us first briefly illustrate the angular moment analysis outlined in Section |TTT| for this 
transparent one-angle distribution. In this case, we have 

g (1 \9) = 3 -(l + cos 2 9), g^ 2 \9) = 3 -sin 2 9. (73) 

Consequently, if we choose 

w {1) (9) = 5 cos 2 9- 1 and w {2) {9) = 2 - 5 cos 2 9, (74) 
the orthogonality relation 

J + ^d{cos9)w^{9)g ( - k \9) = 5 lk (75) 

is satisfied, and we obtain immediately 

l^oC^)! 2 + |-4||(^)| 2 oc ^(5cos 2 ^ - 1) (76) 

\Ai_{t 3 )\ 2 oc ^(2-5cos 2 ^) , (77) 
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where the summation is over all the events in the same time bin as tj. 

In the case of the untagged one-angle distribution, the Amt oscillations proportional 
to the CP-violating weak phase 50 cancel, and the terms ( |76D and (|77|) evolve like 
(|A (0)| 2 + |A||(0)| 2 )e" rL * and |-4j_ (0) | 2 e _rjfff , respectively. A fit (now with only one 
parameter in each time evolution) gives the decay widths Ti and of the CP-even and 
CP-odd B s mass eigenstates, as well as the CP-even and CP-odd rates |A (0)| 2 -|- |A||(0)| 2 
and |v4j_(0)| 2 , respectively. For limited statistics, one may want to use time moments |48| 

roc 

/(«) = / dtt n f(t) . (78) 



The weighting-functions method is thus an alternative to the two-bin method suggested in 
|I|, [£J . Note that we do not need any a priori information about the relative magnitudes 
of CP-even and CP-odd amplitudes. 

In the case of tagged measurements, the integrated decay rates 

r(f)= [ +1 d(cos9)^l (79) 



-i dcos9 
evolve in time for intitially present B s and B s mesons as 

T{t) oc (|A (0)| 2 + |A,|(0)| 2 ) e- FLt + \A ± (0)\ 2 e- THt 

- (\A (0)\ 2 + |A|,(0)| 2 - \A ± (0)\ 2 ) e- Tt sin(Amt) 50 (80) 

and 



T(t) oc (|A (0)| 2 + |A|,(0)| 2 ) e- TLt + |A ± (0)| 2 e- r ** 

+ (|A (0)| 2 + |A,|(0)| 2 - \A ± (0)\ 2 ) e~ Tt sin(Amt) 50, (81) 



respectively, where we have used Tables || and £|. Consequently, the time-dependent CP 
asymmetry arising in the decay B s — > J/ip takes the following form: 

V V ' T(t) + T(t) 

e~ rt sin( Amt) 50 . (82) 





A (0)\ 2 + \ 


4l 


(0)| 2 -|- 


<4±(0)| 2 


(1 


Mo] 


)| 2 + 


^ii (o; 


•I 2 ) 


I e~ T ^ + 


|A ± (0)| 2 e" r ^ 



Using the quantitative estimates collected in Table [I], we obtain 

l^(0)| 2 



|A (0)| 2 + |A||(0)| ; 



0.1... 0.5. (83) 



Although these estimates suffer from large hadronic uncertainties, they indicate that it 
may not be justified to neglect the CP-odd contributions proportional to |^4j_ (0) | 2 in the 
time-dependent CP asymmetry (B3). 
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Observables: 






Mt)\ 2 

A\\(t)\ 2 


i[5(cos 2 #-sin 2 flcos2y?) - 1] 
|[5(cos 2 fl + sin 2 flcos2^) - 1] 
2-5 cos 2 6 


Re(^(t)A|,(t)) 
lm(Al(t)A ± (t)) 
lm(A*(t)A ± (t)) 


sin(2?/>) sin(2y?) 
— | sin(2#) simp 
sin(2?/0 sin(20) cos y2 



Table 5: A set of weighting functions for extracting the observables b^(t) of the decays 
B s -> J/V»(-> and B -> J/V>(-> l + l~)K*(-+ ixK). 

The coefficient of sin( Amt) 5<p in ( |82| ) can be experimentally determined [for instance, 
from the untagged studies outlined above]. Thus the fundamental weak phase 5(j) can 
be cleanly extracted once the Amt oscillations are resolved. Future experiments at the 
Tevatron and the LHC should be able to achieve this goal. Once the Amt oscillations 
are traced, one can alternatively perform a tagged, one-angle, time-dependent study to 
separate the CP-even and CP-odd contributions, from each of which 8(f) can be directly 
extracted. The efficient extraction of the various observables depends on the detector 
configuration, so that other possible variations should be considered. The full angular 
distributions contain, of course, all the available information, and will be determined 
eventually. 

In order to determine 8<p from untagged B s — > J/ip <fi decays, where the Amt oscil- 
lations cancel, the observables corresponding to the interference terms Im (A^(t)A±(t)) 
and Im (AQ(t)A±(t)) must be studied. Valuable information about CP-conserving strong 
phases can also be obtained, thereby sheding light on the hadronization dynamics of 
B s — > J/ip(j) and the issue of "factorization", which predicts trivial strong phases. A set 
of weighting functions applicable to this case is given in Table [5[ 

IV-E The decay B -> J/4)(-> l + l-)K*(^ tcK) 

The angular distribution for B — ► J/ijj{— > l + l~)K*{-^ ttK) takes the same form as 
Eq. (^) if we use the decay angles specified in Eqs. ([65]) and ([66]) with replaced by 
K* and K + replaced by the strange meson. 

Using the same angles for B — > J/i/)(—* l + l~)K (— > irK), we obtain the analogous 
angular distribution to the B s — > J/i/)(—> l + l~)(f)(-^ K + K~) case given in Eq. (|70|). 
The same weighting functions (see Table [5]) can therefore be used to determine the 
corresponding observables in those decays. The comparison of the observables in these 
two modes would give us an idea of the extent of SU(3) breaking. 

If the K*° is observed to decay to the CP eigenstate ir°Ks, the time evolution of the 
corresponding three-angle distributions [Eqs. (|64D and (fTOp l is given in Tables || and 
[7|, respectively Tables |6] and [7] assume that the unmixed amplitudes depend on a 
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Observable 


Time evolution 




Mt)\ 2 
Mt)\ 2 


|A (0)| 2 e- r * [1 + sin(2/3) sin(Amt)] 
U||(0) \ 2 e~ rt [1 + sin(2/3) sin(Amt)] 
\A ± (0)\ 2 e- rt [1 - sin(2/3) sin(Amt)] 


Re(A*(t)A {l (t)) 
lm(Ap)A ± (t)) 
lm(A*(t)A ± (t)) 


A (0) 1 1 Ay (0) | cos(5 2 - 5i)e- r * [1 + sin(2/3) sin(Amt)] 

A||(0) p ± (0)|e" r * [sin(5i) cos(Amt) - cos(2/3) cos(5i) sin(Amt)] 

\Ao{0) \\A±(0) \e~ rt [sm{5 2 ) cos(Amt) - cos(2/3) cos(<5 2 ) sin(Ami)] 



Table 6: Time evolution of the decay Bd — > J/ip{-^ l + l )K*°{— > tt°Ks) of an initially 
(i.e. at t — 0) pure -B^ meson. 



Observable 


Time evolution 




Ao(t)\ 2 
A\\(t)\ 2 
A±(t)\ 2 




A (0)\ 2 e- rt [1 - sin(2/3) sin(Ami)] 
A||(0) 2 e" r * [1 - sin(2/3) sin(Amt)] 
A ± (0)| 2 e~ r * [1 + sin(2/5) sin(Amt)] 


Re(Al(t)A l{ (t)) 
lm(Al(t)A ± (t)) 
1m(Al(t)A x (t)) 




A (0)||A||(0)| cos(5 2 - 5i)e" rt [1 - sin(2/3) sin(Amt)] 
-|A||(0)||A ± (0)|e- rt [sin(5i) cos(Amt) - cos(2/3) cos(<$i) sin(Amt)] 
-|A (0)||A ± (0)|e- rt [sin(5 2 ) cos(Amt) - cos(2/3) cos(5 2 ) sin(Amt)] 



Table 7: Time evolution of the decay B d — > J/ift(—* l + l )K (— > tt .?^) of an initially 
(i.e. at t = 0) pure meson. 



single, unique weak phase, which is justified within the CKM model (see Section II). In 
these tables, V and Am > describe B^-B^ mixing. They are related to each other 
through the mixing parameter Xd = (Am/T) B . In analogy to Eq. QS7D, mixing-induced 
CP violation in B d -> J/V>(-> l + l-)K*°(-^ n°K s ) 
which is given by 



exp(-2i/3) 



49| measures a weak phase /?, 

(84) 



Within the Wolfenstein expansion |T(J , /3 is equal, to a very good approximation, to the 



angle (3 of the "standard" (non-squashed) unitarity triangle j44| . Therefore we have not 
distinguished between j3 and (3 in Tables || and |^. 

Whereas the rates for tagged Bd — > J/ipKs and Bd — > J/ipKs events, which are 
given by 



T[B d (t) -> J/V^s] oc |A(0)|V r ' [1 - sin(2/3) sin(Amt)] 
r[Sd(t) -> J/ipKs] oc |A(0)| V r * [1 + sin(2/3) sin(Amt)] 



(85) 
(86) 



allow only the determination of sm(2/?) and of (Am, T) Bd , an analysis of the tagged 
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three-angle distribution for the decay Bd — > J/^>(— > l + l~)K*°(^ n Q Ks) (and its CP- 
conjugate) yields valuable additional information from the interference terms, as can be 
seen by looking at Tables || and [7]: 

• Ke(AQ(t)A»(t)) provides additional information on cos(^2 — Si). 

• lm(A*^(t)A±(t)), Ihi(Aq(/!;) A±(t)) provide additional information both on sin5i(2) 
and cos<5i(2) and on cos(2/3). The latter quantity plays an important role to resolve 



discrete ambiguities in the determination of the CKM angle (3 [50 



Predictions for these observables are given in Table |l| 

The largest data sample for Bd — > J/ipK*° is, however, not for K*° — > n Ks, but 
for K*° — > ir~K + . The complete angular distributions and time dependences for the 
relevant decay modes are given in Appendix A. For charged B decays, the corresponding 
time and angular distribution is obtained by going to the isospin - related mode and 
setting Am = 0. Experimental studies of these decays are very important, since they 
probe sin(5!( 2 )) and non-factorizable terms through the observables corresponding to the 
left-hand sides of Eqs. (|38|) -(|40"|) |T7j , [19| , [20|[ . The relevant information about 5i and 



<5 2 extracted from these B data samples, "tagged" at the time of decay, can be related 
to B s — > J/ip4> by using SU(3) flavour symmetry of strong interactions, and allows a 
determination of 5<p from the time evolution of even the untagged rate given by Eq. ([TT|) 
||. This approach does not involve the assumption of factorization, just SU(3) symme- 
try arguments. Unfortunately the corresponding S"t/(3)-breaking corrections cannot be 
treated in a quantitative way at present. 

V The angular distribution of the colour-allowed de- 
cays B s D* + D*~ and B -> D* + D* 

The decay of D*^ 1 is predominantly electromagnetic, i.e. D*^ — > Df'j, whereas D* 
decays also strongly to Dir. Therefore the angular distributions of the two decay modes 
discussed in this section are quite different from each other. The first step in the decay 
chain is, however, still of the form P — > V1V2 (as in B s — > J/ipcf)) and consequently 
the terms Aq,A\\,A± retain the same meanings as in Eq. (RBp and the same physical 



significance as in Sect. IV 



V-A The decay B s -> D*+{^ D+-y)D*-(-^ D; 7 ) 



Applying the same convention as in Ref. [51], we define the coordinate system as follows: 



in the rest frame of the decaying B s meson, let the directions of motion of D* + and D*~ 
be z' and z", respectively. In the plane transverse to z' (or z"), choose any direction as 
y' and y". The directions of x' and x" are then specified uniquely via x' = y' x z' and 
x" = y" x z". Thus, x' and x" point in opposite directions. Then (6', if') is the direction 
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w^(9',9", x ) 




A \\ 
A ± 


2 
2 
2 


(-45/92)0(9', 9") + (245/92) sin' 9' sin 2 9" 
(1O/23)C(0',0") - (45/46) sin 2 & sin 2 0" + (25/4) sin 2 9' sin 2 9" cos2x 
(1O/23)C(0', 0") - (45/46) sin 2 9' sin 2 0" - (25/4) sin 2 9' sin 2 0" cos2x 


Re(A* A ll ) 
Im(A|A ± ) 
Im(A*A ± ) 


25y / 2sin0 / sin0" cos0' cos0" cosx 

-(25/4) sin 2 0' sin 2 0" sin2x 
— 25v / 2sin0 / sin0" cos0' cos0" sinx 



Table 8: A set of weighting functions for B s — > D* + (-^ D+-f)D* (— > D77). Here 
C(0', 0") = (1 + cos 2 9')(l + cos 2 0"). 



of Z)+ in the rest frame of D* + in the (x' — y' — z') coordinate system, whereas (0", tp") 
is the direction of D~ in the rest frame of D*~ in the (x" — y" — z") coordinate system. 

Since the choice of directions of y' and y" was completely arbitrary, only the com- 
bination x — <fi' + V 9 " °f V 9 ' an d V 9 " is physical and these two angles will appear in the 
angular distribution only through x- I n terms of the momenta of particles, the angles 
9', 9" and x can be defined as: 

cos9' = p' D + -p D *+ , cos0" = p" D - -p D *- 

sin 9' sin 9" cos x = — cos 0' cos 0" — p'^ • p'' ' D - 

sin 0' sin 9" sin x = (p' D + x p" D - ) • p D *+ . (87) 

The bold-faced quantities are unit three-vectors. The unprimed quantities are measured 
in the rest frame of B s , single-primed quantities in the rest frame of D* + , and double- 
primed quantities in the rest frame of D*~ . 

In terms of these angles, the angular distribution takes the form 

d " T oc JL ( 2 |A | 2 sin 2 9' sin 2 9" 

a cos 9' a cos 9 ax 647T l 

+^|A||| 2 [(l + cos 2 / )(l + cos 2 ,/ ) + sin 2 / sin 2 0" cos2x] 

+ ^|A ± | 2 [(l + cos 2 0')(l + cos 2 0") -sin 2 0' sin 2 9" cos2x] 
— Im (AjjAjJ sin 9' sin 9" sin2x + 2V2 Re (A^A\\) sin 9 sin 9 cos 9 cos 9' cosx 

-2\/2 Im (A* Q A±) sin 9' sin 9" cos 9 cos 9" sin X } , (88) 

where the time dependence of all observables is implicit. It can be read off from Table [| 
The weighting functions are listed in Table [8]. 

The angular distribution for the CP-conjugate process B s — > D* + (— > D+rfD* - (— > 
-D~7) is given by 
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dcosO' dcosO" d\ 6Att 
+ i|^| 2 [(l + cos 2 ^)(l + cos 2 ^")+sm 2 ^ sin 2 #" cos2x] 

+-|Al| 2 [(l + cos 2 00(1 + cos 2 9") - sin 2 9' sin 2 9" cos2 X ] 
— Im (t4|| A±) sin 2 9' sin 2 9" sin2x + 2V2Re (A *A {1 ) sin 9' sin 9" cos#' cos 9" cosx 
- 2 V2 Im (AT^I) 

sin#' sin^" cos#' cos 9" sin X } • (89) 

The time evolution of the various quantities is the same as in Table [|. 

As in the case of B s — > J/i/xj), the decay B s — > D* + D*~ is dominated by a sin- 
gle weak amplitude. Therefore the analysis of the tagged and untagged decays out- 
lined in Sect. |1V-B| and [1V-C] remains valid by replacing (J/ip, 4>,l + ,l~ , K + , K~) — > 
(D* + , D*~ , Df, 7, D~ , 7). Since this process is colour-allowed, factorization is expected 
to hold more strongly. 

V-B The decay B -> D* + (-^ D+-f)D*{-^ Dtt) 

Whereas the decay of the D* + meson is of the form V — > P 7 , which has the same angular 
dependence as V — > Z + /~ for massless leptons, the .D decay belongs to the category 
— > P\Pi- The net angular distribution should therefore have the same form as that for 
P — > V(— > l + l~)V(—> P1P2). The angular distribution is thus given by Eq. (|64|), where 
the definitions of angles are the same as in Eqs. (|65|) and fl66|) with replaced by D , 
Z + replaced by Df, and replaced by the charmed meson arising from the D decay. 
The angular distribution for the CP-conjugate decay B — > D*~(^ D~j)D*(-^> Dtt) is 
given as in Eq. (170). The weighting functions collected in Table [5] can be used to extract 
the corresponding observables from experimental data. 

At this point, a few comments concerning the time evolution of these angular dis- 
tributions are in order. Let us first consider decays of neutral B d mesons. Since here 
the final states are flavour-specific, no interference effects between B d -B d mixing and 
decay processes arise in this case. Consequently, the time evolution of the corresponding 
observables is only governed by the "mixing" of the initial particle, which is either a pure 
B d or B~ d . For B d -> D*+(-^ D+ 7 )D-(-> Dtt) and B~ d -> Dj>y)B* + (-> Dtt), 

the time evolution of a generic observable Q of the angular distribution takes hence the 
form Q(t) = Q(0)e' rt cos 2 (Amt/2), while we have Q(t) = Q(0) e~ rt sin 2 (Amt/2) in 
the case of B d -> D*~(-^ Djj)D*+(-^ Dtt) and B~ d -> £>* + H D+j)D*~(-^ Dtt). The 
time evolutions of the untagged flavour-specific decays and the related B ± decays can 
be obtained straightforwardly from these expressions by setting Am = 0. 
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VI Summary 



The kinematics of B and B s meson decays into two vector-particles, which both con- 
tinue to decay through CP-conserving interactions into two lighter particles, involve three 
independent decay angles. The time evolution of the coefficients of the corresponding an- 
gular distributions contains valuable information about the lifetime and mass differences 
between the B s mass eigenstates Bf and B^, the relative magnitudes and phases of CP- 
odd and CP-even decay amplitudes, and CP-violating effects, including the Wolfenstein 
parameter 77 and the CKM angle f3. The ratios of these coefficients are estimated by 
using various form-factor models. Determinations of these time-dependent coefficients 
will be useful in testing these models and furthermore in determining the extent to which 
factorization or the SU (3) flavour symmetry of strong interactions hold in these decays. 

The observables of the angular distributions can be determined from experimental 
data by an angular-moment analysis in which the data are weighted by judiciously chosen 
weighting functions in order to arrive directly at the observables. At times, this permits 
the extraction of the fundamental CKM parameters. A method applicable to all kinds 
of angular distributions is indicated, where the weighting functions can be determined 
without any a priori knowledge of the values of the coefficients. This method is almost 
as good as the likelihood-fit method for a small number of parameters and is expected to 
give some reliable results even with low statistics where a likelihood fit to a large number 
of parameters is inefficient. 

The B s meson decays B s — > J/ip(/>, D* s + D* s ~ are considered in the light of a possible 
width difference (AT)b s - The observables of their angular distributions can be related 
to those of the decays B — > J/ipK*, D* + D by using the 577(3) flavour symmetry, where 
B stands for B^ or B + . The full angular distributions for all these transitions are given 
explicitly, and the corresponding weighting functions are specified. The time-dependent 
observables in all these decays provide information about the corresponding values of 
Ar and Am. In addition, the decays of B s mesons inform us about the Wolfenstein 
parameter 77, while Bj — > J/ipK*(^ k°Ks) probes the CKM angle (3. Some of the 
quantities related to the B s case can even be extracted from untagged data samples, where 
one does not distinguish between initially present B s or B s mesons. The comparison 
between coefficients of angular distributions of B s and B mesons may give us an idea 
about S77(3)-breaking effects, while the comparison of colour-suppressed (B — > J/ipV) 
and colour- allowed (B — > D^V) modes should help in testing the expectation that 
factorization holds to a greater extent in the latter case. 
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Appendix 



A Angular distributions and time dependences for 
flavour-specific B d — > J/ip(-^ l + l~) K*{-^ K ±f K T ) modes 

(-) 

The angles are defined as in (pq), where the meson is replaced by K*, and the K + 
meson by the strange meson in the final state. In order to parametrize the corresponding 
angular distributions, we use the following combinations of trigonometric functions: 

fx = 2 cos 2 ip (1 — sin 2 9 cos 2 <p) 

fa = sin 2 ip (1 — sin 2 9 sin 2 p) 

fa = sin 2 ip sin 2 9 

fi = sin 2 ip sin 29 sin 93 

fa = (1/V2) sin 2^ sin 2 #sin2</? 

fa = (1/V2) sin 2^ sin 20 cos^. (90) 

Taking into account \Af \ = \ Af \ and using the notation Af = Af(0), where / G {0, ||, _L}, 
we obtain 



d 3 T[B d (t) ■ 


-> J/V>(-> l + l~)K*( 


-> K+7T-)] 9 




dcos9 dp d cos ip 


327T 


x {fx\A \ 


2 + fa\A l{ \ 2 + fa\A± 


| 2 -/ 4 Im(Afi4j> 


d 3 T[B d (t) - 




i^-7T + )] 9 




dcos9 dp d cos ip 


327T 


x {fx\A \ 


2 + fa\A {l \ 2 + fa\A ± 


l 2 + / 4 im (a;a ± )- 


d 3 T[B d (t) ■ 


-> j/v>(-> 


-> 9 




dcos9 dp d cos ip 


327T 


x {fi\A \ 


2 + fa^] 2 + fa\A ± 


| 2 + / 4 Im(^A ± )- 


d 3 r[B d (t) - 


-> J/V>(-> l + l~)K*{ 


^K + 7T-)] 9 




dcos9 dp d cos ip 


327T 



COS 



2 (^) (91) 



cos' e"« (92) 



sin=(^)e-« (93) 



si, 2 (^) (94) 
{/i|A,| 2 + /2|A||| 2 + / 3 |2L| 2 - / 4 Im(XJyl ± ) +hRe(A' A t ) + /elm^x)} . 
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